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Nonparallel Stability of Two-Dimensional Nonuniformly 
Heated Boundary-Layer Flows 

A. H. Nayfeh and N. M. £1-Hady 

Engineering Science and Mechanics Department, 

Virginia Polytechnic Institute and State University, 

Blacksburg, Virginia 24061 

An analysis is presented for the linear stability of water bound- 
ary-layer flows over nonuniformly heated flat plates. Included in the 
analysis are disturbances due to velocity, pressure, temperature, den- 
sity, and transport properties as well as variations of the liquid 
properties with temperature. The method of multiple scales is used to 
account for the nonparallelism of the mean flow. In contrast with 
previous analyses, the nonsimilarity of the mean flow is taken into 
account. No analysis agrees, even qualitatively, with the experimental 
data when similar profiles are used. However, both the parallel and 
nonparallel results qualitatively agree with the experimental results of 
Strazisar and Reshotko when nonsimilar profiles are used. 



I. Introduction 


Linke^ was the first to investigate experimentally the effect of 

heat transfer on transition. He measured the drag on a vertical heated 

plate that is placed in a horizontal airstream. He found that heating 

the plate causes its drag to increase considerably. He concluded from 

this observation that heating the plate is destabilizing. Liepmann and 
2 

Fila fully confirmed the destabilizing effect of heating in air bound- 
ary layers. They performed measurements on a vertical heated plate in a 
horizontal airstream. They found that the critical Reynolds number 
decreases with wall heating. The destabilizing effect of heating in an 
incompressible air boundary layer is due to increasing the air viscosity 
next to the wall, thereby producing an inflected velocity profile. On 
the other hand, cooling yields a fuller velocity profile and hence a 
more stable flow. 

Since heating decreases the viscosity of water, the above measure- 
ments and arguiP.ents suggest that heating the surface of a body in a 

water stream is stabilizing. This has been confirmed by the analysis of 
3 4 

Wazzan, et al ’ . Their analysis is for a parallel flow and is based on 
the disturbance-vorticity equation only; that is, it does not include 
the energy equation and hence the temperature fluctuations. However, 
their analysis includes the effects of the mean-temperature distribution 
on the viscosity of the fluid. With these assumptions, Wazzan, et al^’^ 
obtained a fourth-order modified Orr-Sommerfeld problem. Their results 
show that the critical Reynolds number increases as the wall heating 
increases, reaches a maximum, and then decreases. Lowell and Reshotko^ 
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reformulated the parallel stability problem and included temperature as 

well as vorticity fluctuations. They ended up with a sixth-order rather 

than a fourth-order system. They found that the solutions of the fourth' 

and sixth-order systems are close for all wall temperature ratios over 

the normal liquid range of water. 

The stabilizing effect of small wall -temperature ratios in water 

was confirmed experimentally by Strazisar, et al® and Parker^. Parker 

found that the transition Reynolds number for water flowing in a tube 

can be increased from 10 x 10^ to 42 x 10® by using a 7'’C wall overheat. 

Strazisar, et al® conducted an experiment for the case of uniform wall 

overheat. Their results show that, as the wall heating increases, the 

critical Reynolds number increases, the growth rates decrease, and the 

range of frequencies undergoing amplification decreases. All of these 

3-5 

results qualitatively agree with the parallel stability results . To 
compare quantitatively with the experimental results, El-Hady and 

O 

Nayfeh used the method of multiple scales to develop a nonparallel 
stability theory for heated water boundary layers. The nonparallel 
results are in good agreement with the experimental data. 

Since the flow over the portion of the body upstream of the criti- 
cal Reynolds number is stable, no stabilization is needed on that por- 
tion, and one would need heating only on the portions downstream of the 

critical Reynolds number. This suggests the use of nonuniform rather 

g 

than uniform wall heating. This led Strazisar and Reshotko to examine 
experimentally the effect of nonuniform wall heating. They conducted 
experiments with two types of wall heating. The first is a step change 
in temperatures and the second is a power-law temperature variation of 
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the form T - T = Ax", where T is the temperature, x is the distance in 
w c 

the streamwsie direction, n is a constant, and the subscripts w and e 

denote conditions at the wall and the edge of the boundary layer, respec' 

tively. In their power-law case, they kept ~ fixed while 

they changed the exponent n. They made all their measurements at 

which corresponds to a displacement- thickness Reynolds number of about 

800. Their results show that decreasing n is stabilizing; that is, the 

case n < 0 results in lower growth rates than the case n = 0 (uniform 

case), which in turn results in lower growth rates than the case n = 1. 

These results could not be explained, even qualitatively, by using the 
9 10 

parallel analyses ’ . This led to the speculation that an appropriate 
nonparallel theory may be needed to explain these results. However, 
applying the nonparallel theory with a similar mean flow, we were also 
unable to explain, even qualitatively, these results as shown in Fig. 1. 

Looking closely at the aforementioned parallel and nonparallel 
calculations, one finds that all of them employ self-similar boundary-, 
layer profiles. For a uniform wall temperature or for a power-law 
temperature distribution in a fluid having constant properties, the flow 
is self-similar. However, for a fluid with variable properties, the 
flow is not self-similar if the wall temperature is not uniform. In 
fact, the mean-flow measurements of Strazisar and Reshotko show varia- 
tions of the mean flow from the similar solution. Therefore, the pur- 
pose of the present paper is to examine the parallel and nonparallel 
stability of nonsimilar water boundary layers over nonuniformly heated 
flat plates. 



II. Problem Formulation and Method of Solution 


The present study is concerned with the two-dimensional, nonparallel 
stability of two-dimensional, viscous, heat conducting liquid boundary 
layers to small amplitude disturbances. The analysis takes into account 
variations in the fluid properties but neglects buoyancy and the dissi- 
pation energy. All fluid properties are assumed to be known functions 
of the temperature alone. Dimensionless quantities are introduced by 
using a suitable reference length L* and the frees tream values as refer- 
ence quantities, where the asterisk denotes dimensional quantities. 

To study the linear stability of a mean boundary-layer flow, we 
superpose a small time-dependent disturbance on each mean-flow, ther- 
modynamic, and transport quantity. Thus, we let 

q(x,y,t) = + q(x,y,t), (1) 

where Q^{x,y) is a mean steady quantity and q(x,y,t) is an unsteady 
disturbance quantity. Here, q stands for the streamwise and transverse 
velocity components u and v, the temperature T, the pressure p, the 
density p, the specific heat Cp, the viscosity p, and the thermal 

conductivity k. Substituting Eq. (1) into the Navier-Stokes and energy 
equations, subtracting the mean quantities, and linearizing the result- 
ing equations in the q's, we obtain the following disturbance equations: 

If + (PsU +pUs) +|y (PsV + pVs) = 0, (2) 
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an an 3U 9U 3U^ 9U 

Ps^lt ^ ^ ^ ^ ''s fc V'' + p(^s 

= - 1^ ^ 1 [ir §7 ^ ^ 


+ ^ [y (|ii + 1^) +w(:r-2. + .^)]] 

3y ^^s'3y 3x' ^'3y 3x '-'■ 


3U. 3V. 


^s^3t ■*■ ^s 3x “ 3x ■*■ '^s 3y 3y ^ ■*■ 3x^ * '^s 3y^^ 

= - |£ + 

K |^3X "^S'3y 3X 

3V 


9v . .. 


3v . .. ^''s, 


3V. 


(3) 

3V. 


3y 


a an a„ 3U, 3V^ 


3U. 




j- 


(4) 

3T. 


^ ^ ^ ^s IJ ^ 3T ^ ''s ly^ ^(Ps ^ 3x 


iL 91 + . !Is-, ^ , 31 

s 3y ■‘ RP»*gCp 3x 'S 3x ^ 3x ■' 3y '^s 3y 


1 

+ V v^] = 

■ c ^ \/ J 




P. U. K. Cp = 


dOs dp^ d<^ '^Ppj 


dig ’ dTj ’ dT^ » dT^ 


T . 


(5) 

( 6 ) 


Here, c is the liquid specific heat at constant pressure, R = 

P$ 

p*U*L*/p* is the Reynolds number, and Pr = c* p*/<* is the freestream 
e e e e p^ e e 

Prandtl number. Moreover, 


r = I (£ + 2), m = I (i - 1). f = ^ (1 + 2£), X*| U-1), 

(7) 
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where 1 Is the ratio of the second to the first viscosity coefficients 
(£ = 0 is the Stokes assumption). 

The problem is completed by the specification of the boundary 
conditions; they are 


at y = 0, 

(8) 

as y -*• «>. 

(9) 


We restrict our analysis to mean flows that are slightly nonparallel; 
that is, the transverse velocity component is small compared with respect 
to the streamwise velocity component. This condition demands all mean- 
fiow variables to be weak functions of the streamwise position. These 
assumptions are expressed mathematically by writing the mean-flow vari- 
ables in the form 

Ug = UjCxi.y), Vj = eVj{x,,y), = P5(xi), T^ = Tj(xi.y), 

(10) 

where Xj = cx with e being a small dimensionless parameter character- 
izing the nonparallelism of the mean flow. In what follows, we drop the 
carrot from V^. 

To determine an approximate solution to Eqs. (2)-(10), we use the 
method of multiple scales^^ and seek a first-order expansion for the 
disturbance variables u, v, p, and T in the form of a traveling harmonic 
wave; that is, we expand each disturbance flow quantity in the form 

q(xi,y,t) = [qo(xi.y) + eqi(x,,y) + ...]exp(ie), (11) 

where 

If = a.i(xi). If = - 0 ). (12) 
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For the case of spatial stability, Oo is the complex wavenumber for the 
quasi-parallel flow problem and oj is the disturbance angular frequency, 
which is taken to be real. 

Substituting Eqs. (11) and (12) into Eqs. (2)-(10), transforming 
the time and spatial derivatives from t and x to 0 and Xi, and equating 
the coefficients of and e on both sides, we obtain problems describ- 
ing the qo and qi flow quantities. These problems are referred to as 
the zeroth- and first-order problems and they are solved in the next two 
sections. 
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III. The Zeroth-Order Problem 


Substituting Eqs. (11) and (12) into Eqs. (2)-(l0) and equating the 
coefficients of e" on both sides, we obtain the following problem: 

Li(uo,Vo,Po,To) = iooEOjUa + (U^ - ^)Po] + (pjVo) = 0, (13) 


3U. 


L2(uo,Vo.Po.To) = [iP5ao(Uj - ^) + ^ Psao]Uo + (Pg 


i S I 0 o t % S\ I S dUg 

R 3y ■ R 3y 'dT. 3y ^ " R 3y 3y 


if 


3Vn 


■ R ^s“° 3y 


1 ^^s 3Tn 


3^Uo 


D Pr ^ 


R dT^ 3y 3y R ^s 3y 


L3(uo,Vo,Po,To) = [iPgCXo(Ug - “) + PgCto]Vo “ ^ “o 


(14) 


3u, 


- R ao 


T if 3Ug 

dig 3y ■ R 3y 


3Vi 
R 3y 3y 



(15) 


L4(uo,Vo,Po,To) = [iPgaodJj " o7^ 


1 


RPr_c 


e Pc 


3^< 


^ 2 

3y ■ RPr„c„ 3y 3y " RPr^c 


1 .. ^ _ 


e Pc 


e Pc 


<0 “ P* (15) 


The boundary conditions are: 

Uo = Vo = To = 0 at y = 0, 
Uo, Vo, To 0 as y -> «. 


(17) 

(18) 


Equations (13)-(18) constitute an eigenvalue problem, which is 
solved numerically. It is convenient to express it as a set of six 
first-order equations by introducing the new variables z defined by 

on 
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N 

0 

tl 

C 

0 

Zq2 

_ aufl. 

~ 3y ’ 

Z03 

= Vo, 

Zo4 - Pd« 

0 

N 

= To. 

2o6 

II 


Then, Eqs. (13)-(18) can be rewritten in the compact form 


9Z .. 

g 



clL 

8y 

- 1 a. .z . = 0 
j=l 

for i = 1,2.., 6, 

(20) 

Zci = 

Zo3 - Zos “ 0 

at y = 0, 

(21) 

2oi » 

Zo3« Zo5 ■* 0 

as y 

(22) 


where the a^j are the elements of a 6 x 6 variable-coefficient matrix. 
The nineteen nonzero elements of this matrix are listed in Appendix I. 

To set up the numerical solution, we first replace the boundary 
conditions (22) by a new set at y = y„ where y„ is a convenient location 
outside the boundary layer. Outside the boundary layer, the mean flow 
is independent of y and the coefficients a. . are constants. Hence, the 

* J 

general solution of Eqs. (20) can be expressed in the form 


z 


oi 


6 

I A. .c.exp(A.y) for i = 1,2, 
j=l 


.6. y = yg 


(23) 


where the X. are the eigenvalues of the matrix [a..], the A., are the 
J 

elements of the corresponding eigenvector matrix, and the c. are arbi- 

J 

trary constants. The real parts of three of the X. are negative, while 

•J 

the real parts of the remaining X. are positive. Let us order these 

eigenvalues so that the real parts of Xi,X2, and X3 are negative. Then, 

the boundary condition (22) demands that C4, C5, and Cs are zero. To 

set up this condition for the numerical procedure, we first solve Eqs. 

(23) for the c-exp(X.y) and obtain 
J J 
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(24) 


6 

c.exp(A.y) = I b -.2 . for j = 1,2,... 
J J ^-*1 *J 01 

where the matrix [b. .] is the inverse of [A. . 

I J I J 

= 0 in Eq. (24) leads to 
6 


I = 0 for j = 4,5, and 6 at y 


. 6 , 

]. Setting C4=Cs=Cs 

" V 


(25) 


Using Eqs. (25) as the boundary condition at y = yg and guessing a 

value for oo, we integrate Eqs. (20) from y = yg to y = 0 by using the 

12 

comouter program developed by Scott and Watts that employs a Gram- 
Schmidt orthonormal ization procedure, and then we attempt to satisfy the 
boundary conditions (21). If the guessed value for qo is the correct 
eigenvalue, the three boundary conditions will be satisfied. In general, 
the guessed value is not the correct value and the boundary conditions 
at the wall are not satisfied. A Newton-Raphson procedure is used to 
update the value of oo and the integration is repeated until the wall- 
boundary conditions are satisfied to within a prescribed accuracy. This 
leads to a value for «o and the eigenfunctions are recovered using the 
stored solution vectors. They can be expressed in the form 

z . = A(xi)r.(xi,y) for i = 1,2,.. .,6, (26) 

0 < * 

where A is still an undetermined function at this level of approxima- 
tion. ^t is determined by imposing the solvability condition at the 
next level of approximation. 
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IV. The First-Order Problem 


With the solution of the zeroth-order problem given by Eq. (26), 
the first-order problem becomes 


3Z . 6 

v-^-- - y a. .z . = G. 
3y ij iJ 1 

^+D^A for i = 1,2 6, 

(27) 

Zll = Zl3 = Zi5 = 0 

at y = 0 , 

(28) 

Zll, Zi3, Zi5 -*• 0 

as y , 

(29) 


where the and are known functions of the oo. and the mean-flow 
quantities. They are defined in Appendix II. 

Since the homogeneous parts of Eqs. (27)-(29) are the same as Eqs. 
( 20 )-( 22 ) and since the latter have a nontrivial solution, the inhomo- 
geneous Eqs. (27)-(29) have a solution if, and only if, a solvability 
condition is satisfied. In this case, the solvability condition demands 
the inhomogeneities to be orthogonal to every solution of the adjoint 
homogeneous problem; that is, 

G Jfl 

i, Sr * 

where the W. (xi,y) are the solutions of the adjoint homogeneous problem 
corresponding to the eigenvalue oo. Thus, they are the solutions of 

9W. e 

j aj.Wj = 0 for i = 1,2,...,6, (31) 

Wz = W 4 = Ws = 0 at y = 0, (32) 

Wa, W 4 , Wg -*■ 0 as y -*• “>, (33) 

Substituting for the G. and D. from Appendix II into Eq. (30), we 
obtain the following equation for the evolution of the amplitude A: 
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(34) 


A dxi " 


where 


■ - y J, 

•'o /) '* . 

The solution of Eq. (34) can be written as 

A = Aoexp[ie ^ai(xj)dx], 
where Aq is a constant of integration. 


(35) 


(36) 


To determine aj(xt), we need to evaluate doo/dxi and the 3;;^/3xi. 
To accomplish this, we differentiate Eqs. (20)-(22) with respect to Xj 
and obtain 


I7 - J, It * f""- < = 

) at y = 0, (38) 

! as y (39) 

The initial conditions for the computational procedures are chosen to 
exclude any multiple of the homogeneous solutions. The S^- and G^. are 
known functions of oo. and the mean-flow quantities and their deri- 
vatives; they are given by 




3Ci 

= ^ = 

3Cs 

3xi 

3Xi 

3xi 

3Ci 

3Cs 


3x 1 

’ 3Xi ’ 

3xi 


6 3a . . 

s. = y c. — ^ 

i 3xi 

J- 1 


and G. 


3a,, 


"i "" -L S' ’ ~ 1 , 2 ,..., 6 (40) 

3 ^ 

Using the solvability condition of Eqs. (37)-(39), we find that 
^ G,W,dy], 

•'n ^0 


dg 

dx 


(41) 
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Therefore, to the first approximation 

zoi = AoCi(xi,y)exp[i j [cto + eai)dx - iojt] + 0(e), (42) 

where the are related to the disturbance variables by Eq. (19) and 
the constant Aq is determined from the initial conditions. It is clear 
from Eq. (42) that, in addition to the dependence of the eigensolutions 
on Xi, the eigenvalue Qq is modified by caj. The present solution re- 
duces to those obtained by Nayfeh, et al^^ and Saric and Nayfeh^^ for 

I 

the case of nonheat conducting flows. 
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V. The Mean Flow 


For flows whose thermodynamic and transport properties are func- 
tions of temperature, the two-dimensional boundary-layer equations for 
a zero-pressure gradient and in boundary- layer coordinates are 


I 7 (pu) + ^ (pv) = 0 , (43) 

^ 3y 

9u . — 3u 3 / 3U\ 

puc + pvc ^ ^ (k . (45) 

P P 3y 3y 3y 


The temperature dependence of p and u couples the momentum and 
energy equations. Note that buoyancy and viscous dissipation effects 
are neglected. 

We introduce the Levy-Lees transformation^^ 


dC = PgUgPgdx , (46) 

dn = pUg(2C)’'/'dy . (47) 

Then, the derivatives with respect to x and y are transformed according 
to 


3 _ II 3 , 3r| 3 

3x ^ Pe^e^e 3? 3x 3n ’ 


37 ' (20 


Substituting Eqs. (48) and (49) into Eqs. (43)-(45) yields 


(48) 

(49) 
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(50) 




where 


U = u/Ug, H = T/Tg , 


V = 


2C 


p U U 
^e e^e 


3n 

ox 2(5) /!■ 




C = py/pgUg, Pr = CpP/< . 


The boundary conditions are 


U(C.O) = 0. V(C.O) = 0, H(C,0) = H^(C). 
U(c,n) 1 and H(?,n) -*■ 1 as n Hg. 


(51) 

(52) 

(53) 

(54) 

(55) 

(56) 

(57) 


Equations (50)-(52), (56), and (57) are numerically integrated using a 
step by step procedure in the streamwise direciton. A three-point 
implicit finite-difference technique is used to reduce the energy and 
momentum equations and the boundary conditions to a set of simultaneous 
tridiagonal equations. These equations are linearized and then solved 
by using the algorithm of Thomas. Then, the continuity equation Is 
numerically integrated by using the trapezoidal rule. The method of 

solution closely parallels those of Fliigge-Lotz and Blottner^®, Davis 

17 18 

and Fliigge-Lotz , and Harris . 
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VI. Analytical Results and Comparison with Experiments 

g 

Strazisar and Reshotko performed their experiments in a water 
tunnel whose test section was 394 mm long, 229 mm wide, and 152 mm high. 
The freestream turbulence was 0.1 - 0.2% for U* < 3.4 m/s. They mea- 
sured the boundary- layer characteristics on a flat. plate that was 348 
mm long and 16 mm thick and spanned the test section. The plate was 
fitted with a rounded leading edge (0.79 mm radius) located 10.8 mm 
below the top of the test section. 

Disturbances were artificially introduced in the boundary layer by 
using a vibrating ribbon that is stretched across the plate surface 95.3 
mm behind the leading edge. The amplitudes of the generated distur- 
bances were measured at five stations spaced 6.4 mm apart between x = 

127 mm and X = 152.4 mm. They traversed the boundary layer in the 
normal direction and recorded the peak amplitude. Then, they determined 
the growth rates at x = 139.7 irm by using a polynomial curve fit of the 
peak amplitude data. 

The plate heating was provided by 11 electric heaters distributed 
along the plate. The wall temperature was monitored by using 11 ther- 
mistors imbedded in the surface of the plate at its centerline. How- 
ever, because of the large temperature gradients involved, the ther- 
mistors did not accurately yield the plate temperature. Consequently, 
they had to determine the wall temperature from boundary- layer profiles 
measured with a hot-film anemometer operating as a resistance thermo- 
meter. Due to equipment limitations, the wall temperature could not be 
monitored or maintained near the leading edge. Since the thermal bound- 
ary is very thin near the leading edge, measurement of temperature 
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profiles using the hot- film anemometer were impractical in that region 
and the first measurement of the wall temperature was provided by a 
thermistor imbedded 30.5 mm from the leading edge. 

In the case of the power-law distribution T - T » Ax", Strazisar 

we 

and Reshotko held the temperature difference fixed at x^.g^ , while they 
varied n and x as shown in Fig. 2. They presented growth data at x^g^ 

= 139.7 mm only (corresponding to R = 475). Their results show that 
decreasing n is stabilizing. However, Fig. 1 sliows that decreasing n is 
destabilizing in both the parallel and nonparallel calculations when a 
similar mean flow is used. Hence, using self-similar mean profiles 
cannot predict the experimental data. 

Since the upstream wall temperature distribuiton is essential for 
calculating nonsimilar boundary layers we are unable to compare quanti- 
tatively the analytical results with the data of Strazisar and Roshotko 
for the case of power-law distributions. Figures 3-5 show the variation 
of the parallel and nonparallel growth rates with frequency (defined as 
F = ii)*v*/U*^) calculated for the power-law distributions shown in Fig. 2 
for AT = 1.67®, 2.78°, and 4.44°C at x^g^ = 139.7 mm. In each figure, 
we show the results for n = -0.5, 0, and 1 as well as the results for 
the unheated case. The nonparallel growth rates do not include the 
distortion effect of the mode shape. Including this distortion modifies 
quantitatively but not qualitatively the results. Both parallel and 
nonparallel theories predict that decreasing the exponent n results in a 
stabilizing effect at x^g^ , in qualitative agreement with the experi- 

g 

mental results . 

The stabilizing effect produced by decreasing the exponent n can be 
explained as follows. As n decreases. Fig. 2 shows that AT increases 
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at all locations upstream of . But increasing AT results in a 

fuller velocity profile and hence a more stable flow. Therefore, the 
stabilizing effect produced at is a cumulative of all upstream 

stabilizing effects. However, as n decreases, AT decreases downstream 
of , resulting in less fuller velocity profiles. Therefore, at 
some location downstream of x^^^ , a distribution with a larger exponent 
will be more stabilizing as shown in Fig, 6. Thus the neutral stability 
curves are not nested and conclusions regarding stabilizing and desta- 
bilizing effects away from depend on the Reynolds number. The 
integration of the growth rates yields the amplification factor, which 
seems to be the best indicator of stability. Figure 7 shows the varia- 
tion of the integrated growth rates of Fig. 6 with Reynolds number. 

Also shown is the variation of the maximum amplification factor with the 
exponent n. It appears that the maximum amplification factor for n = 

- 0.5 and n=0 are nearly the same, while it is higher for n = 1. This 
result holds for other frequencies as indicated by the nesting of the 
growth-rate curves (Fig. 4) as functions of frequency. 
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VII. Conclusion 


We analyze the linear nonparallel stability of two-dimensional 
liquid boundary layers on a flat plate for the case uf nonuniform wall 
heating. Stability calculations using a self-similar mean flow cannot 
predict, even qualitatively, the experimental results of Strazisar and 
Reshotko for power-law temperature distributions. However, by using 
nonsimilar mean flows, both parallel and nonparallel results are in 
qualitative agreement with the experiments. The stabilizing and desta- 
bilizing effects at x^.^^ depend on the temperature distribution. 
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APPENDIX III 


The variation of the thermodynamic and transport properties with 
temperature is given by Lowell and Reshotko^ 


P* = 


1 - 


(T* - 3.9863)2(1* + 288.94U1 
508929.2(T* +68.12963) 


+ 0.011445 exp(- 


j.4. 3 

T* 


), 


p* in gm/m£, T* in "C. 


= 1.37023(7* - 2 0) + 8.36 x 10‘^(T - Pnl* 
U* ' 109 + T* 

P* in Cp T* in “C 

K* = - 0.901090 + 0.1001982T* - 1.873892 x 10“ 

+ 1.039570 X 10' V’ , 

K* in megawatts cm“^K'\ T* in °K. 


cj = 2.13974 - 9.68137 x 10" V + 2.68536 x 10" V* 
- 2.42139 X 10" 

cj in cal gm"\"\ T* in "K . 
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Figure Captions 


Figure 1. Variation of the spatial amplification rate with frequency 
for power-law temperature distributions for 
AT = 2.78°C at using similar mean-flow profiles. 

Figure 2. Power-law temperature distributions. 

Figure 3. Variation of the spatial amplification rate with frequency 

unheated case and the power-law tempera- 
ture distributions of Fig. 2 for AT = 1.67°C at x , using 

I C I • 

nonsimilar mean-flow profiles. 

Figure 4. Variation of the spatial amplification rate with frequency 
^ref. unheated case and the power-law temperature 

distributions of Fig. 2 for AT = 2.78°C at x^^^ using non- 
similar mean-flow profiles. 

Figure 5. Variation of the spatial amplification rate with frequency 
^ref, unheated case and the power-law temperature 

distributions of Fig. 2 for AT = 4.44»C at x^^^ using non- 
similar mean-flow profiles. 

Figure 6. Variation of the spatial amplification rate with streamwise 
position for the unheated case and the power-law temperature 
distributions of Fig. 2 for AT = 2.78°C at x^^^ using non- 
similar mean-flow profiles. 

Figure 7. Variation of the amplification factor with streamwise posi- 
tion for the unheated case and the power-law temperature 
distributions of Fig. 2 for AT = 2.78°C at x^^^ using non- 
similar mean-flow profiles. 
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Figure 1. Variation of the spatial amplification rate with frequency 
at for power-law temperature distributions for 

AT = 2.78°C at using similar mean-flow profiles. 
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Figure 3. Variation of the spatial amplification rate with frequency 


at for the unheated case and the power-law tempera- 
ture distributions of Fig. 2 for AT = 1.67°C at using 
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Fiqure 4. Variation of the spatial amplification rate with frequency 
at for the unheated case and the power-law temperature 

distributions of Fig. 2 for AT = 2.78°C at using non- 

sini]ar mean- flow profiles. 



Figure 5. Variation of the spatial amplification rate with frequency 


at X r for the unheated case and the power-law temperature 
ref. 

distributions of Fig. 2 for AT = 4.44°C at using non- 


similar mean-flov/ orofiles. 
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Variation of the spatial amplification rate with stre^im • 
position for the unheated case and tho nr, i streamwise 
distributions of Fig. ^ for AT = ? temperature 

similar mean-flow profiles. 
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Figure 7. Variation of the amplification factor with streamwise posi- 
tion for the unheated case and the power-law temperature 
distributions of Fig, 2 for M = 2.78'’C at using non- 

similar mean-flow profiles. 
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